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Abstract 

Bound muons (also pions, kaons, etc) increase the fission barrier and produce some stabilizing 
effects for highly charged nuclei. If the binding energy of the muon exceeds mc 2 , it becomes stable. 
The Is state of a muon inside an exotic nucleus with atomic number A = 5Z/2 and such large 
charge Z that the Is energy E is in the range < E < —mc 2 is considered. 
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As well known, all atoms with nuclear charge Z > 100 have very short lifetimes and 
atoms with Z > 118 have not as yet been found. However, the behavior of orbitals in the 
field of a highly charged "nucleus" remains an interesting theoretical problem. There are 
various motivations for the study of such systems. Firstly, in collisions of heavy nuclei (e.g. 
uranium with charge Zj = 92) a compound nucleus with the charge Z = 2Zt is formed 
for a short time. Electron-positron pair creation due to "vacuum breakdown" is expected 
in this situation, provided the total compound nucleus charge Z exceeds a critical charge 
Z c fa 170 (see e.g. review [lj and the book 2|]). An experimental investigation of this 
phenomenon has been done in GSI, with new experiments in preparation jjj. Then, there 
are speculations about exotic forms of matter such as "strange matter" where a nucleus 
contains strange quarks in addition to "normal" u and d quarks Q]; there may be highly 
charged small primordial black holes j^] and the associated black hole "atoms" j^]. 

In this note, we would like to consider the purely theoretical problem of a charged particle 
in the superstrong Coulomb field. The electron problem is actually very well studied (see 
e.g. l(]). However, the Coulomb problem for a heavier particle has some new and interesting 
features. Consider, for example, a muon. Normally, it decays into electron and neutrinos. 
However, the bound muon is stable if the binding energy exceeds mc 2 so the total relativistic 
energy E < 0. In this case the bound muon has a significant stabilizing effect on the 
"nucleus", it contributes ~ mc 2 ~ 100 MeV to the nuclear fission barrier. In the case 
of bosonic particles like pions or kaons the effect can be even more significant since many 
particles may be placed in Is state and the critical charge for spin S = and S = 1 is much 
smaller than for S — 1/2 (Z c =68 for a very small nucleus, i.e. a factor of two smaller than 
for Dirac particles). Moreover, the strong interaction may produce an additional attraction 
and reduce the critical charge. 

In our numerical calculations, we assume a standard nuclear charge density with radius 
R = 1.15A 1 / 3 fm, and atomic number A = 5Z/2. Let us start our consideration from the 
conventional electron case. The electron density is localized mostly outside the nucleus even 
as the nuclear charge approaches the critical value Z c fa 170; thus, there is no collapse of the 
electron wave function inside the nucleus. In Fig. [H we plot the large and small components 
of Is wave function and the radial density for charge 2 = 166, which is close to the critical 
charge. The electron energy in this case is E = —0.796 mc 2 (the energy for the critical charge 
is E r = —mc 2 ). Thus, the electron wave function for the finite nucleus in the main area 
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FIG. 1: Radial wave functions P(r) and Q(r) and the radial charge density Airr 2 p(r) of a Is 
electron (E = — 0.796 mc 2 ) in the field of a nucleus of charge Z = 166 are shown. 

is not very different from the point-nucleus case. Consequently, the critical charge for the 
finite nucleus Z c « 170 is relatively close to the case of a very small nucleus where Z c = 138. 

Now consider the muon case. The non-relativistic muonic atom is m M /m e = 200 times 
smaller than the usual electronic atom. In the relativistic case E < 0, the muon wave 
function for the Is state is localized inside the nucleus where the electrostatic potential has 
the oscillator shape 

V (r) = — 



R \2 2B? 

The critical charge for this oscillator potential is very different from the Coulomb potential. 
The energy reaches E = for a nuclear charge Z w 840 and the critical value E = —mc 2 
for Z c w 2000. 

In the non-relativistic case, the radius of Is state is inversely proportional to the nuclear 
charge Z. Surprisingly, in the relativistic case E < 0, the muon wave function is practically 
independent of Z (see Fig. [21 upper panel). The explanation is as follows: In the main 
area, the electrostatic potential V(r) is approximately constant. A constant potential V 



shifts the energy (by V), however, it does not change the wave function since the Dirac (and 
Schrodinger) equations contain the combination E — V which does not change. To illustrate 
this explanation, we plot E — mc 2 — V(r) for Z = 840 (E « 0) and for Z = 2000 (E « — mc 2 ) 
in the lower panel of Fig. [2] The corresponding wave functions practically coincide (see the 
upper panel). In Fig. [31 we present the radial densities of Is, 2s, 3s and 4s states for critical 
charge Z c « 2000. All of these states are localized inside the nucleus. 

Note that the "critical charges" for muons are so large that there is hardly any chance 
of actual stabilization of a "normal" nuclear matter by muons. Therefore, in the muon case 
this problem is only an exercise in the relativistic quantum mechanics (critical charges in 
the oscillator-like potential of the finite nucleus). However, this solution may give us some 
hints as to what one should expect in more interesting cases of pions and kaons. We do not 
consider these cases in the present note since they require more sophisticated calculations 
including the strong interaction. 
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FIG. 2: Upper panel: The Is radial density for a muon (energy slightly above — mc 2 ) in the field of a 
nucleus with charge Z = 2000 and radius R = 19.2 fm (red curve) is seen to be essentially identical 
to the Is radial density of a muon (energy 0) in the field of a nucleus with Z = 840, R = 14.6 fm 
(black curve). Lower panel: The functions E — mc? — V(r) for the E = —mc? case (red curve) and 
E = case (black curve) are essentially identical in the region where the densities are nonvanishing, 
explaining the fact that the densities agree in the two cases. 
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FIG. 3: Upper panel: Radial densities for the four lowest ns states of a muon in the field of a 
nucleus with charge Z = 2000 and radius R = 19.2 fm. Lower panel: Energies of the lowest four 
ns states are shown in comparison with the potential energy V(r). 
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